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I. INTRODUCTION 

The recent work by Komargodski and Schwimmer [1, 2] on the a theorem has attracted much 
interest. The a theorem is thought to be the generalization to four dimensions of Zamolodchikov's 
c theorem [3] which establishes the irreversibility of the renormalization group (RG) flow in two 
dimensions. Several extensions of the work of Zamolodchikov to four dimensions have been 
proposed. In particular, Cardy pointed to the function a given by the integral of the trace of the 
energy-momentum tensor over the four-sphere [4]. He then conjectured that a has properties 
similar to its two dimensional cousin c. Cardy's conjecture was tested within the framework of 
perturbation theory in [5-7]. It was discovered that the a function is not monotonically decreasing 
along the RG flow, and an alternative function, denoted here by a, was devised to be monotonically 
decreasing in perturbation theory. A nice property of a is that it coincides with Cardy's a function 
at fixed points, where the jS functions vanish. 

Komargodski and Schwimmer suggested a non-perturbative proof of the a theorem by intro- 
ducing a dilatonic field [1, 2] whose scattering amplitude is related to the a function. In fact, the 
analyticity of the amplitude supports the monotonicity of the a function (not a) along the RG flow 1 . 
Establishing the existence of a function that is monotonic along the RG flow can lead to relevant 
constraints on the dynamics of a given gauge theory. In supersymmetry the existing relation 
between a and the R-charge [10, 11] found use in the a-maximization [12]. Using holographic 
methods, the a theorem could be tested and verified in the context of supersymmetry [13, 14]. The 
generalization, using holography, to arbitrary space-time dimensions was explored in [15]. More 
recently, the relation between scale and conformal invariance in 4D was elucidated in the context 
of the a theorem [16, 17]. 

The goal of this work is to extend the perturbative analyses of the a theorem for nonsuper- 
symmetric gauge theories with fermions and gauge singlet scalars, interacting via Yukawas, to 
the maximum known order in perturbation theory. This allows us to investigate the details of the 
a theorem, particularly for non-standard fixed point structures. 

The paper is structured as follows. In Section II we introduce the essential tools and discuss 
the generic expression of the a function relevant for its determination to three loops in the gauge 
coupling. Concrete examples are presented in Section III where we explicitly evaluate a at the 
Banks-Zaks (BZ) fixed point (FP) for a gauge-Yukawa theory to the leading and the next-to- 
leading order in perturbation theory. We then discover, for a certain non asymptotically free gauge 

1 The relation between the trace anomaly and the amplitude of the three point function of the energy-momentum 
tensor was already pointed out in [8, 9]. 
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theory, the appearance of a perturbative UV fixed point (UVFP), and at the next-to-leading order 
also the appearance of another infrared (IR) fixed point. By changing the number of matter fields 
we observe the merging and disappearance of both fixed points. This rich structure of fixed points 
constitutes an interesting playground for elucidating the properties of the a function in gauge 
theories. We conclude in Section IV and provide a number of technical details in Appendix A. 

II. THE a THEOREM BEYOND THE LEADING ORDER 

In four dimensions and for a general quantum field theory the vacuum expectation value of 
the trace of the energy-momentum tensor for a locally flat metric g f(V reads 

(r;;) = cW 2 (^ F )-fl£ 4 (^v) + -- - , (i) 

where a and c are real coefficients, E^(g^ v ) the Euler density and W(g^ v ) the Weyl tensor. The dots 
represent contributions coming from operators that can be constructed out of the fields defining 
the theory. Their contribution is proportional to the f> functions of their couplings. The coefficient 
a is the one used in Cardy's conjecture, and for a free field theory it is [18] 

1/11 \ 
90(8tz) 2 V 2 ; / w 

where n s , n t and n v are respectively the number of real scalars, Weyl fermions and gauge bosons. 

The change of a along the RG flow is directly related to the underlying dynamics of the theory 
via the jS functions. This can be shown by exploiting the abelian nature of the trace anomaly 
which leads to the Weyl consistency conditions in much the same manner as the well known 
Wess-Zumino consistency conditions [19]. As discussed in the introduction and following the 
work of Jack and Osborn [5, 6], rather than using a one uses the function a related to it by 

a = a + Wfii , (3) 

where W,- is a one-form which depends on the couplings of the theory. The Weyl consistency 
conditions 2 imply for a 

dfi = -Xifij + (diWj - djWtfj , (4) 



2 It is worth mentioning that the Weyl consistency conditions used above assume that the trace of the energy- 
momentum tensor vanishes when all the jS functions are zero simultaneously. Exceptions are known to exist [20] and 
in this case one modifies the consistency conditions [16, 17] in order to build a. 
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where Xij can be viewed as a metric in the space of couplings. The positivity of the metric x is 
established in perturbation theory, and therefore in this regime the function a is monotonic along 
the RG flow 

Pj^ =Pidia = -XijPiPj<0. (5) 

The irreversibility of the RG flow has been conjectured to be valid beyond perturbation theory. 

In the next subsections, capitalizing on the work of Jack and Osborn [6], we will construct a for 
gauge theories featuring fermionic matter interacting with gauge singlet scalar fields via Yukawa 
interactions to the highest known order in perturbation theory. 



A. Gauge- Yukawa Theories 

We consider the following Lagrangian skeleton 

£ = ~^, lv P lv + iWyD^i + \d^ a d^a ~ V?V } <t>a + h.C.) - ^A abcd ^b^d , (6) 

where we dropped the gauge indices for F^ v and the Weyl fermions W,. The fermions transform 
according to a given representation R of the underlying gauge group. The real scalar s (p are singlets 
with respect to the gauge group. The indices run over the number of matter fields. To exemplify 
our results we will consider gauge theories for which the Yukawa and quartic interactions depend 
each on a single parameter as follows 

Vlj = V ^ii ' ^abcd = A Tabcd , (7) 

where the T are tensors specified in a given theory. There are therefore three couplings in our 
setup: gauge g, Yukawa y and quartic A. We define 

(4n) z (47z) z (4n) z 

The generic structure of the associated jl functions {§ a = da/dln^i) reads 3 

fia g = -2a 2 g [fe + hocg + bydy + b 2 a 2 g + b 3 a g a y + ha 2 ,] , (9) 
p ay = 2a y \c\a.y + c 2 a g + c 3 a g a y + c^a 2 + c 5 a 2 + C(,a y a A + c 7 a^] , (10) 
Pa A = d 1 a 2 A + d 2 a A a y + d 3 a 2 . (11) 

3 The factors of 2 in the definition of the gauge and Yukawa jS functions follow from the definitions (8). One has for 
example %Ja t = 2fi g /g. 
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The expansion to three loops in the gauge coupling, two loops in the Yukawa and one in the 
quartic coupling leads to a consistent expression for a to order a 3 . If the scalars were charged 
under the gauge group, terms proportional to a g a\ would appear. 

Having established the generic form of the jS functions, we move to determining the metric x 
and one-form W. They can be found by examining the relevant Feynman diagrams which enter 
the computation of the trace anomaly, as shown in Appendix A. We find 



X 



%f (l + Aa g + Bia z g + B 2 a g a y ) B 

Bo ^ (l + B 3 a y + B,a g ) 

Xaa ) 



(12) 



The coefficient Xgg enters at the one-loop order, A and Xyy at two loops, while xaa and the B;'s 
appear only at three loops. Similarly, the one-form W takes the form 



W g - — (Dq + D\a g + C\a z + C2«gay) , 



W y = D 2 + C 3 a v + C 4 a g , (13) 
W A = D 3 a A . 

The general structure of x confirms that it is sufficient for all our purposes to consider the Yukawa 
jS function (10) to two-loop order and the quartic one (11) to one-loop only. 
The leading coefficients Xgg> Xyy an d Xaa are [6] 

d(G) 1 11 \ 1 / 1 T \ ( . 

X SS ~ 128ti 2 ' Xyy ~ 1287Z 2 \3 V V/' XAA ~ 1287I 2 1 72 bcd bcd j ' ' 

where we used the tensors T defined in eq. (7) and d(G) denotes the dimension of the adjoint 
representation G of the underlying gauge group, i.e. the number of gluons. A is given by [6] 

A = 17C 2 (G)-yN R T(R), (15) 

where C 2 (G) is the quadratic Casimir of the adjoint, Nr the number of Weyl fermion in the 
representation R and T(R) is the trace normalization satisfying T(R)5 a b = Tr(R a R},), R a being the 
generators of the fermions under the gauge group. With these coefficients, the metric x is positive 
definite near the origin of the coupling constant space. It is however clear that in the absence of a 
theorem, the positivity of x away from the origin is not guaranteed. The remaining coefficients of 
X and W are yet to be determined but, as we shall show, they are not needed to determine a at the 
fixed points to the order investigated here. 
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B. Power of the consistency relations and the a-function 

The system of first order differential equations in (4) allows to derive the following conditions 
relating the different coefficients of the $ functions as well as x and W, 



\/ = -XyyCl , XyyCe = Xxxh , ^XyyCj = X\\&1 , 
2-Xggh + Xyy (B4C1 + B 3 c 2 + c 3 ) = 2 (B - C 2 + C 4 ) Ci , (16) 
(B 2 fc + Ab y + b 3 ) + 2* w (B 4 c 2 + c 4 ) = (B - C 2 + C 4 ) c 2 + 2 (B + C 2 - C 4 ) £> • 

The equations in the first line above can be used to either test or determine some of the higher 
order coefficients of the jS functions since we know the metric factors. The remaining equations 
can be used in a similar fashion. However, given that the B, coefficients have not been explicitly 
computed we use the knowledge of the jS functions to deduce, for example, B3 and B 4 assuming 
that c 2 does not vanish. 

For a to cubic order in the couplings, and using the consistency relations above, we have 

a = fl free + fl (1) + fl (2) + fl (3) + . . . , (17) 

where flf ree is the free-field theory value (2), and the one, two and three-loops coefficients are 

= -2xggha g , (18) 
a (1) = -x g g (h + Ab ) a Z g ~ 2xg g b y a g a y + XyyCia 2 y , (19) 



« (3) - Xgg 



- (b 2 + Mi) a 3 g + (b 3 + Ab v ) a\a y + 2b 4 a g a y + ~ ^4& 4 ~^( h 3+ Ab y fj a y 
I (cs - ^c 3 + ^ J c 4 J a y + c 6 a 2 y a A + 2c 7 a y a\ + ^Xaa«i«a ( 20 ) 



+ Xyy 

fia g r 1 3 \ ^a./ Bq - C 2 + C 4 



+ 



ly 4c 2 
where we defined 



/(«?) = *«(f«J + - 1 (If 4) ♦ (If 4 ■ PI) 

Remarkably, the unknown coefficients B, and Q appear only in the last two terms of aS 3 ', where 
they are multiplied by p 1 functions and hence vanish at fixed points 4 . This was also observed to 
occur in super symmetric theories [21]. 



4 We have taken the liberty of adding higher order terms in order to rewrite the coefficients as jS functions. These terms 
are irrelevant when computing a between perturbative fixed points. 
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It is instructive to calculate a to the second order using (3) and recalling that the leading 
coefficients entering the one-form W are Do = Xggr D\ = \ J&L Xgg an d D 2 = \Xyy [6]. We find the 
simple expression 

a = fl free + Xgg (ha 2 g + b y a g a y ) + O (of) / ( 22 ) 

where, remarkably, the term linear in «g, the one quadratic in ay as well as the term linear in A 
canceled out. Because the signs of b\ and by depend on the gauge theory, a is not generally a 
monotonically decreasing function along the perturbative RG flow. 



C. a at fixed points 



We can now move to the study of the fixed points and determine the variation of a between two 
of them. A convenient way to search for the zeros of the system of jS functions (9)-(ll) is to first 
solve analytically for = which permits to relate a a to a y , then set to zero jS a further relating 
cty to a g , so that finally we can search for the zeros of the following effective jS function in a g 



?f g = -2a 2 g [bo + bfa g + bfal], 



(23) 



where 



with 



bf = b X ~ -by , 

1 C\ y 



c 2, 



z C\ \C2/ C\ C\ \C\) 



„eff 



eff d 3 d 2 ( d 2 

^ =C5 -d 1 C7 -2dA C6 -ck C7 



1+ Jl- 



(24) 
(25) 



(26) 



At this order in perturbation theory, there can be at most two perturbative fixed points for each 
sign in c^ ff , if both bo and b^ are tuned to be small. An example of this is provided in the following 
section. Using eq. (24) and (25), the difference in the function a — or equivalently a — between 
the UV and IR fixed points can then be written as 



Aflperturbative = (^ ~ ^perturbative = -2X g g\b ~ «f ) + r (bf + Ab ) ({aff - {aff) 



+ 



(bf + Abf + Bbo^a^f-iaff) 



(27) 
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where a^ v and a g R denote the values of the gauge coupling at the UV and IR fixed point respec- 
tively, and we defined 

BsBl _£2( B2+ M. (28) 



c i \ Xgg) 

This expression reduces to the case of a gauge theory without Yukawa interactions by replacing 
bf and bf with b\ and bi. 

Inspecting the effective jS function, there is a perturbative fixed point for small bo, which reads 

«f = -^+0(»j). (29) 

For positive bo and negative bf , this is the usual Banks-Zaks IR fixed point. The situation in 
which the BZ fixed point is of UV nature is equally possible. This occurs by reverting the signs 
of both bo and bf . The trivial fixed point at the origin will be in the first case an UV fixed point 
and in the second an IR one. The finite change in a between the UV and IR fixed points can be 
computed either way, and one obtains 

MBZ = ■ (30) 

Here the sign reflects the sign of bo, and Aa is then positive for any physical fixed point. However, 
A5 can formally become negative when the value of the coupling a g at the fixed point is on the 
unphysical negative axis. 

To three loop order in the effective jS function, one can have the two following physical zeros 

bf 



a BZ = 



°T L , gpg \ bf/ 4P p r ^ 



4Z> Q bf\ 



For small values of bo, the solution with negative sign corresponds to the usual BZ fixed point, 
with the following corrections 

a BZ = _^( 1+ b ^f +0 ( h Z\)_ (32) 
g bf{ (bf) 2 y0, j 

This expression holds provided bo/(bf) 2 is small. We shall see below that there are cases where 
this is not true. Using (32), we can compute the three loop corrections to the variation of a, 



M™ = *XggMl-iAbf-2bf) : /! '< 



(33) 



bf{ 1 2 3(bf) 2 j ' 

We now turn our attention to the second zero a^r. The first observation is that for a generic value 
of bf this fixed point occurs at a non-perturbative value of the coupling. This is what happens 
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in general for gauge theories with fermionic matter in a given irreducible representation of the 
gauge group [22]. However, for gauge theories with Yukawa interactions and/or multiple matter 
representations, the possibility that both bo and b^ f are small exists. An explicit example is provided 
below. Furthermore, when b~ - (b^) 2 / (4bo) the two fixed points merge. This phenomenon can 
happen within the range of perturbation theory. At the merger, one has ^merger = -2bo /b^, which, 
when plugged into Eq. (27), gives 

M BZ | =*Xg£\- (34) 

I merger 'a*3^erf \ / 

The virtues of this expression will be studied in more detail elsewhere [23]. 

Having in our hands the explicit tools, we can explore the a theorem for gauge theories with 
interesting fixed point structures. 

III. A CONCRETE EXAMPLE 

We consider a SU(N C ) gauge theory with Nf fundamental Dirac fermions Q = I adjoint 

Weyl fermions A, and a gauge singlet complex scalar H that transforms in the bifundamental 
representation of the SU(Nf)L X SU(Nf)R global symmetry of the theory. For the benefit of the 
reader, the field content and the quantum symmetries of the theory with t — 1 are summarized in 
Table I. The Lagrangian of the theory is 

1 



£ = Tr 



2 -F^F fIV + mx + QilDQ + d^dm + y H QHQ 



ill (Tr [H r H]) z - u 2 Tr(H H) . (35) 



Tr is the trace over both color and flavor indices. is the usual covariant derivative. Throughout 



TABLE I. Field content of the example. The first three fields are Weyl spinors in the (|,0) representation 
of the Lorentz group. H is a complex scalar and G„ is the gauge field. U(1)af is the extra Anomaly Free 
symmetry arising due to the presence of A. 



Fields 


[SU(N C )] 


SU(N f ) L 


SU(N f ) R 


U(l) v 


U(1) AF 


A 
1 


Adj 

□ 
□ 


1 

□ 
1 


1 
1 

□ 




Nf-N c 

N c 
Nf-Nc 

N c 


1 

N c 
Nf 

N c 
N f 


H 


1 


□ 


□ 





2N C 
Nf 


G, 


Adj 


1 


1 
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this section we will work with the rescaled couplings which enable a finite Veneziano limit of the 
theory with i fixed. That is, we let N c ,Nf — > oo while keeping x = Nf/N c fixed. The appropriately 
rescaled couplings are 

g 2 N c f H N c u,Nj UlNf 

fl * = (^r flH -(4^- Zl -(4^r Z2 = (4^- (36) 

This model was introduced in [24] to investigate near-conformal dynamics, at the one loop level, 
and its impact on the spectrum of the theory with special attention to the dilaton properties. The 
model was further investigated at the two loop level in [25]. To compute a, following the previous 
section, we need to determine the three loop contribution to the gauge f> function. 
Using [26-31] we find in the Veneziano limit 



fi ag = - - a 2 g [ll - 21 - 2x + (34 - 16£ - 13x) a g + 3x 2 a H 



81x 2 3x 2 (7 + 6x) 2 2857 + 112x 2 - x(1709 - 257/) -1976£ + U5£ 2 2 

+ — a 8 a H 2 U H + is a g 



\, -s r ,n 20(x + ^) - 203 , Q x(x + 12) 2 „ 2 
2(x + l)a H - 6a g + (8x + 5)a g a H + — a g - 8xz 2 a H — — -a^ + 4z^ 



(37) 



(38) 



p Z2 =2(2z 2 a H + 4z 2 2 -xa 2 H ) . (39) 

Here one can see that the double trace coupling z\ does not participate in the running of the 
remaining couplings. In addition, using (14) and (15) the metric coefficients for this theory can be 
found: 

N 2 N j N } 10 

^ = ^2- ^ = 3^?' XAA = 3~2ht' A = 17 ~Y {x + £) - (40) 

One can check that the expressions above satisfy the consistency relations given in the first line 
of (16), and therefore it constitutes an independent check of the correctness of the jl functions. We 
now turn to the FP analysis of the model which will reveal an interesting perturbative structure. 



A. Leading order analysis: Banks-Zaks fixed point 

In order to see a physical BZ fixed point, the one-loop coefficient of the gauge fi function has to 
be small and the signs of bo and have to be opposite. Therefore, our first task is to find a region 
in the parameter space of the model where the physical BZ fixed point exists. We use Eq. (24) 

b = ^(ll - 2{£ + x)) , bf = ^(34 - 16£ - 13x + . (41) 
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0.04 
0.03 
0.02 
0.01 
0.00 
-0.01 
-0.02 
-0.03 



\ Yukawa 
.No Yukawa \ 



4.0 4.2 4.4 4.6 4.8 5.0 

X 

Fig. l.a Aa between the Gaussian and BZ fixed 
points normalised to Xgg at leading order for the 
I = 1 case. The red (dashed blue) line corresponds to 
the model without (with) Yukawa interactions. In 
both cases the physical BZ FP is an IR one. 



A a 



0.03 
0.02 
0.01 
0.00 
-0.01 
-0.02 
-0.03 



■ ^^^^foYukawa^ 


Yukawa / 

/' 
y 


t =U 








/ 



















5.0 



5.2 



Fig. l.b Aa between the Gaussian and BZ fixed 
points normalised to Xgg at leading order for the 
i = case. In the absence (presence) of Yukawa 
interactions, the physical BZ fixed pont is an IR (UV) 
one. The color code is the same as on the left panel. 



From the asymptotic freedom (AF) boundary condition bo = we obtain that x = (11 - 1€)j1. 
Substituting this value of x into b^, we have 



1AF 2 



25 



3(11 - ief 

4^-26 



(42) 



where the last term comes from the Yukawa interactions. We immediately notice that for I* « 0.37 
the b^ F vanishes. Below, we will consider the cases I - 1 for which this coefficient is negative 
and t = for which it is positive. In the first case we have a standard IR BZ fixed point, and 
in the second we obtain a new UV BZ fixed point. It is worth noticing that, in the absence of 
Yukawa interactions, b^ F in (42) is always negative and therefore the physical BZ FP can only be 
the standard IR fixed point. 



1. 1=1 case 



In this case there exists a perturbative IR fixed point regardless of whether we consider the 
presence of Yukawa interactions. We show in Fig. l.a the leading order result for the change in a 
between the Gaussian (trivial) FP and the BZ IR one at leading order, both in the presence (blue 
line) and absence (red line) of Yukawa interactions. Both curves cross zero at x* = 9/2, when 
asymptotic freedom is lost. For x > 9/2, where bo < 0, there is an unphysical BZ UV fixed point 
yielding a negative Aa. The Yukawa interactions in (42) imply \b^ f \ < \b\\, which leads to a larger 
A5 in the case of the gauge theory with scalars. 
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X 



Fig. 2.a The next-to-leading order physical FP 
structure for the 6 = 1 case with Yukawa and quartic 
interactions. 



A a 




x 



Fig. 2.b AS normalised to Xgg for the t = 1 case. The red 
and dashed blue lines are leading order results from 
Fig. l.a while the dotted black and green lines are the 
next-to-leading order corrections. 



2. I = case 



We now turn to the I = case where, rather than having a BZ IR fixed point, the theory develops 
a UVFP when asymptotic freedom is lost, i.e. for bo < 0. Of course, this is possible only because 
of the presence of the Yukawa interactions. In Fig. l.b, the leading order result for the change in 
a between the Gaussian and BZ fixed points without Yukawa interactions is shown in red, and 
the one with Yukawa interactions in blue. Both curves cross zero for x* = 11/2 when asymptotic 
freedom is lost. 



B. Next-to-leading order analysis: Fixed point merger 

At the next perturbative order we deal with the full system of Eqs. (37)-(39) and from now on, 
we concentrate only on the physical fixed points. 



1. £ = 1 case 



We start again with the I — 1 theory and in Fig. 2.a we display the FPs structure for the model 
with Yukawa and quartic interactions. We notice that at x « 3.25 the FP value of the gauge coupling 
vanishes. However, this happens in the region beyond applicability of perturbation theory since 
the two remaining coupling constants are large. In Fig. 2.b we plot the change in the fl-function for 
the next-to-leading order BZ IR fixed point and compare it with the corresponding leading order 
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Fig. 3.a The next-to-leading order physical FP 
structure for the I = case with Yukawa and quartic 
interactions. The vertical dashed-dotted green line 
represents the point where asymptotic freedom is 
lost. 



0.004 




Fig. 3.b AS normalised to Xgg for the I = case. The red 
and dashed blue lines are leading order results from 
Fig. l.b while the dotted black and green lines are 
next-to-leading order corrections. 



results from Fig. l.a. As a general feature, we notice that the next-to-leading order corrections 
reduce the value of Aa in the perturbative regime. It is clear from the plots that for the theory with 
Yukawa interactions perturbation theory breaks down, when moving away from the critical value 
x* — 4.5, earlier with respect to the theory without Yukawa interactions. 



2. I = case 



We now turn to the I - theory where, as discussed above, there is no BZ IR fixed point. 
However, in the asymptotically free regime there is a new physical IR fixed point emerging at 
the next-to-leading order. This non-BZ IR fixed point is present also when asymptotic freedom 
is lost, i.e. for x > 5.5. In this region there is also a BZ UV fixed point that was discussed in 
subsection III A 2. The complete FP structure is shown in Fig. 3. a. The UV and IR fixed points 
merge around x « 5.6. In Fig. 3.b we plot the change in the a-function at next-to-leading order 
together with the corresponding leading order results from Fig. l.b. We notice that Aa becomes 
negative just before the merger which is incompatible with the a theorem. We interpret this effect 
as the breakdown of the perturbative expansion since the FP values of the couplings at the merger 
are quite large, as can be seen from Fig. 3. a. 

So far, all our calculations of Aa were for the flow connecting the trivial FP at the origin of the 
coupling constant space with the BZ one. However, it is relevant also to determine Aa for the 
branch connecting the two non-trivial FPs. In the theory with i — and x > 5.5 this is the RG flow 
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5.52 5.54 5.56 5.58 5.60 



x 

Fig. 4.a A3 normalised to Xgg for the RG flow 
between the BZ UV fixed point and the non-BZ IR 
fixed point. 




5.1501 5.1502 5.1503 5.1504 



x 

Fig. 4.b A5 normalised to Xgg for the i = 0.35 case. The 
magenta curve corresponds to the flow between the BZ 
UV and the non-BZ IR fixed points. The green curve is 
the result between the Gaussian IR fixed point and the 
UV BZ one. 



between the BZ UV fixed point and the non-BZ IR one. We display the change in the fl-function 
in Fig. 4. a. Of course, at the merger Aa vanishes. 



3. t— 0.35 case: The perturbative merger 

Using both i and x as continuous parameters it is formally possible to study the merging 
phenomenon within the perturbative regime. This happens around £ » 0.37. Therefore we 
provide an example with t - 0.35. The change in the fl-function for the two RG flows are shown 
in Fig. 4.b. Since for this value of i perturbation theory holds, we observe a positive and well 
behaved Afl all the way to the merger. 



IV. CONCLUSIONS 



We studied the a theorem to the maximum known order in perturbation theory for nonsuper- 
symmetric gauge theories with fermionic matter as well as gauge singlet scalar fields, where the 
latter two interact via Yukawa interactions. The computation involves three loops in the gauge 
p 1 function, two loops in the Yukawa interactions and one loop in the quartic interactions. To this 
order we have first determined the general expression for the change in the function a between 
two fixed points, and then specialised it to the case of QCD with extra adjoint Weyl fermions and 
elementary mesons interacting via Yukawa terms. We employed the Veneziano limit and deter- 
mined the three loop jS functions of the theory. This has allowed us to test the a theorem beyond the 
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lowest order in perturbation theory. We discovered that the model posses an interesting structure 
of the fixed points, depending on the number of adjoint fermionic matter fields, among which 
we highlight the presence of a Banks-Zaks UV fixed point as well as the occurrence of a merging 
phenomenon for which the a theorem had not been tested before. 
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Appendix A: Diagrammatic form of the trace anomaly 

The structure of the metric (12) and one-form (13) entering the trace anomaly can be determined 
by looking at vacuum polarisation diagrams containing in addition to the usual vertices of the 
quantum field theory the following counterterms 




ft 5 o iZ£ „o o a P"S o P"y o a 

g a? p S a V2 SPg Pa g Pg a m Py fl i/2 pA pa A 

The terms entering the metric (12) are the ones proportional to two powers of the jS functions while 
the ones with one power of jS, fix the one form W (13). All the vacuum polarisation diagrams up 
to three-loop order as well as the form of their contribution are shown in Table II. Note that more 
diagrams would be present if the scalar field were charged under the gauge group. 
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